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We study an ensemble of random matrices (the Rosenzweig-Porter model) which, in contrast to 
the standard Gaussian ensemble, is not invariant under changes of basis. We show that a rather 
complete understanding of its level correlations can be obtained within the standard framework 
of diagrammatic perturbation theory. The structure of the perturbation expansion allows for an 
interpretation of the level structure on simple physical grounds, an aspect that is missing in the 
exact analysis (T Guhr, Phys. Rev. Lett. 76, 2258 (1996), T. Guhr and A. Miiller-Groeling, 
:ond-mat/9702113|) . 

PACS number: 05.45.+b Theory and models of chaotic systems 



I. INTRODUCTION 



Random matrix ensembles were introduced into 
physics by Wigner, Dyson and others Q as phenomeno- 
logical models of complex quantum systems. Such en- 
sembles are defined as to obey certain symmetries and 
constraints but are otherwise "as random as possible" . 
For instance, the Gaussian unitary ensemble (GUE) con- 
sists of all N x N Hermitian matrices H , the only con- 
straint being that, on the average, tr H 2 is a given con- 
stant. This leads to a probability density in the matrix 
space, Vq ~ exp (— tr iJ 2 ) , which is invariant under uni- 
tary transformations. 

Recently there has been some interest in various gener- 
alizations of the GUE and its orthogonal and symplectic 
counterparts 0] . One possible generalization amounts to 
breaking the U(N) symmetry of the GUE by introducing 
a parameter /i into the probability density function 



P(H)dH = Af exp {—F{H))dH , 



(1) 



N 



f(h) = y. h *+ 2 ( 1 +^Y, t( Re H ^ 2 + ( Im ^) 2 ] 



Here Hij, with i < j, designate the independent ma- 
trix elements of a N x N Hermitian matrix, dH = 

N 

\~[ dHu Yl d(Re H^) d(Im Hij) is the volume element 

i—l i<j 

in the matrix space and M is a normalization constant. 
For /i = the expression in the curly brackets is equal 
to (— tr H 2 ), so that the GUE is recovered. The pa- 
rameter fj, breaks the U(N) symmetry and introduces a 
preferential basis. When fi — > oo, for fixed N , all matri- 
ces become diagonal in that basis. The ensemble, thus, 
exhibits a crossover from the Wigner-Dyson statistics of 
the standard random matrix theory (/z = 0) to the Pois- 
son statistics of uncorrelated levels (/i = oo). Such an 
ensemble (for real symmetric matrices) was introduced 
by Rosenzweig and Porter H in their studies of complex 



atomic spectra, and more recently appeared in the field 
of quantum chaos and localization 0- 

We shall be interested in the behavior of this ensem- 
ble in the N — > oo limit. In this limit, significant de- 
viations from the GUE-behavior can only occur if /i in- 
creases with N sufficiently fast. The local level statistics 
is controlled Q - M by the parameter n/N 2 . Only when 
this parameter approaches infinity, the Wigner-Dyson 
statistics becomes completely obliterated and the Pois- 
son limit of uncorrelated levels is reached. In the opposite 
case, i.e. when (fi/N 2 ) — * 0, an arbitrary large sequence 
of levels will obey the Wigner-Dyson statistics of the 
GUE. The "critical" case corresponds to /i = cN 2 , with 
c = const. The situation resembles the one which occurs 
in disordered electronic systems where, in the thermody- 
namic limit, three distinct types of statistics correspond- 
ing to insulator, metal and a mobility-edge system [^]|| 
exist. 

In the present paper we shall take a closer look at the 
eigenvalue statistics, with an emphasis on the "critical" 
case n = cN 2 . We will show that a rather complete 
picture emerges already from diagrammatic perturbation 
calculations, along the lines of Refs. [|l0(ll|. In this case 
the two-point correlation function R(s) (smoothed out 
over few level spacings) differs substantially from both, 
the GUE and the Poisson correlation functions. Here s 
denotes the energy difference in units of the average level 
spacing. For small c, R(s) is approximately given by its 
Wigner-Dyson value, — l/(27r 2 s 2 ), as long as s <C 
p"2[ . For larger s, however, R(s) changes sign, reaches a 
maximum and eventually decreases as l/(7rcs 2 ). 



II. DIAGRAMMATIC ANALYSIS 



To begin with, let us introduce a definition of the 
Rosenzweig-Porter (RP) model which is equivalent to 
Eq. (|l|) but more convenient for diagrammatic compu- 
tation. We define 



1 



H = H Q + V , (H )i 



o 



(2) 



where Ei are independent real random numbers with 
Gaussian distribution 



Vie) 



(3) 



The matrix elements Vij of the hermitean matrix V are 
independent complex random numbers with Gaussian 
distributed real and imaginary parts. The distribution 
is determined by 



(V^) = 0, (\Vij\ 2 ) 



1 



2(1 + M) 



(4) 



It is easy to see that the probability density function of 
the thus defined Hamiltonian is just Eq. 

The diagrammatic analysis amounts to a locator ex- 
pansion of the full single particle Green's function 



H) 



E±i5 S\0 : 



(5) 



with respect to the off-diagonal V. The unperturbed 
(bare) Green's function Go — G is called the locator. 

We consider the density of states p(E) = tr 5(E — H), its 
average value 



u(E) = ((p(E)) v ) £ 
and its correlation function 



R(E,E') = {(p(E)p(E')) 



VI e 



v{E)v{E') 



(6) 



(7) 



Here 



((. . .) ) stands for averaging with respect to 



the off-diagonal elements Vij (the diagonal elements Si). 

We concentrate on energy separations u> = E' — E for 
which both energies are close to the middle of the band, 
i.e. close to E = where v{E) is maximal. In this re- 
gion, v{E) is approximately constant, v{E) = v{E') up 
to 0(1/N) relative corrections. The density-density cor- 
relator R(s) = R(E,E')/v 2 will then be a function of 
the dimensionless level separation s = u>/A where the 
average level spacing is A = l/v{E — 0). 

We next analyze the spectral correlation function R(s) 
in the regime s > 1 where perturbative methods are ap- 
plicable. To begin with, we decompose R(s) according 
to R(s) = Ri(s) + R 2 (s) where 



Ri(s) = ^{ ME + w)) v (p(E))v) e - 
-{{p(E + u)) v ) e {{p(E)) v ) 

R2(s) = -^{((p(E + u;)p(E)) v - 
-{p{E + uj)) v {p{E)) v 



(8) 



(9) 



Note that the decomposition R = R\ + R2 is exact. The 
physical significance of the two functions R1/2 will be 



discussed below. Here we merely note that R\ measures 
correlations remaining in the GUE-averaged density of 
states whereas R2 focuses on the GUE-correlations as 
such. 

Representing the density of states in terms of the single 
particle Green's function 

p(E) = -Tr^Imtr G + (E), 

and making use of the fact that correlations (of any type) 
between products of purely retarded/advanced Green's 
functions vanish for N — > 00: (G +n ) = (G + ) n , we obtain 

Rr = ^Re (tr (G+(E + u)) v tr (G~(E)), 



R-2 



2tt 2 

A 2 

2tt 2 



Re ( (tr G + {E + u)trG-(E) 



V.r 



(10) 



where (. . .) c denotes the connected average, (XY) C = 
(XY) — (X) (Y). Before turning to the actual calculation 
of these functions let us make a few methodological re- 
marks and introduce some building blocks that will be 
of importance throughout. The whole approach will be 
based on a perturbative expansion of the Green's func- 
tions in powers of V. It is instructive to visualize the 
structure of the expansion scheme graphically. To this 
end we introduce the notation of Fig. 1(a) where i and 
j represent matrix indices (which will not be indicated 
explicitly unless necessary). As a first step of our per- 
turbative analysis (cf. the definition of the correlation 
function Ri above) we have to calculate the I^-average 
of the Green's function G. In a diagrammatic language, 
the expansion of the Green's function can be visualized 
as is shown in Fig. 1(b). The subsequent diagrammatic 
analysis of this equation is simplified drastically by two 
observations that hold to leading order in TV -1 . 

i) Contributions with 'crossed GUE-lines' (see 
Fig. 1(e)) are negligible pj. 

ii) Diagram segments which are separated from each 
other by vertices Vij are statistically independent 
with respect to the average over the on-site distri- 
bution functions p{e). 

The second statement is based on the fact that the in- 
dices i and j in a diagram like the one shown in Fig. 1(f) 
are eventually summed over independently of each other. 
(All contributions where one of the summations is con- 
strained will be of higher order in TV -1 .) On the other 
hand, the variables £i at different sites are statistically 
independent. Put together, these two facts imply that ii) 
holds to leading order in N^ 1 . 

The diagrammatic expression for the ^-averaged 
Green's function GiSij = {Gij) v , subject to rule i), is 
displayed in Fig. 1(c). The statement ii) implies that to 
leading order in AT -1 , the self energy part T, ± (as shown 
in Fig. 1(d)) can be replaced by the Q averaged one. We 
thus obtain 



2 



N 



In order to solve this equation self consistently, we have 
to compute the energy average of Q . Anticipating that a) 
the self energy will be largely imaginary ~ =p«r and 
b) T -C 1, we obtain 



1 



r l/2 



cfe- 



J5 - £ ± if 



(E - e + iT)(E - e - iT) 



The second equality is based on the assumption that the 
energy argument £ < 1 is close to the middle of the 
band. As a consequence the real part of the integral is 
of O(E) and negligible in comparison with the imagi- 
nary part 0(1). This justifies the assumption a) above. 
Collecting everything so far we obtain the V^-averaged 
Green's function 



Qt 



i 



E - si ± iT 



r = 



2(1 + 0)' 



(11) 



(Note that assumption b) above holds for all \x ~ N x , x > 
1, i.e. Eq.([Tl|) indeed represents the self consistent solu- 
tion of the Dyson equation in Fig. 1(c).) We next insert 
this expression into the defining equation of the correla- 
tion function R\ and obtain 



Ri 



A 2 N 
~2^ r 



Re 



1 



r l/2 



de- 



1 



E- 



iT E 



iT 



A 2 N 2T 



r 3/2 



U3* 



4r 2 ' 



Noting that the level spacing A = it 1 / 2 /N , we arrive at 
the final result 



Ri(s) 



l 



l 



A 

2T ' 



(12) 



for the first of the above introduced correlation functions. 
We postpone the discussion of this equation until the 
complementary correlation function i?2 has been calcu- 
lated. In principle one might compute Ri via a straight- 
forward perturbative expansion of the Green's function 
G. However, experience gained from the analysis of sim- 
ilar correlation functions |TT| has shown that it is ad- 
vantageous to represent the Green's functions according 
to 

G ± (E) = 8e \n(E ± — H) 

prior to the perturbative expansion. In this way we are 
led to consider 



A 2 9 



E'E 



E<=E+u 



tr \n(E'+ - H) tr ln(£T - H) 



V,c 



Expanding the logarithms in powers of V and applying 
the non-crossing rule i) we obtain 



A 2 

R2 = 2^ 9e '' e 



Re ^-(S n (£,£;')> s , (13) 



E'=E+lu 1 — ' n 
n=2 



where the diagrammatic representation of S n (E,E') is 
shown in Fig. 2. There the segments on the outer (inner) 
ring correspond to the Green's function Q + (E') (Q~(E)) 
and the two rings are connected by n T^-lines. (Note 
that a n — 1 contribution is absent because the poten- 
tial V is off-diagonal in the site indices.) The rule ii) 
implies that each segment of the 'wheel' above can be 
averaged individually over the on-site energies £j. As 
a result, the d factorizes, S n = 7™, 7 = 



JV/(2(1 + fj))(G + (E')g-(E)) E and £ r 



— ln(l — 7) — 7. Computing the energy average (cf. the 
computation of the correlation function R\ above) 



(g + (E')g-(E)) e = 



2v 1 ' 2 i 



E> - E + 2iT ' 



and collecting constants we obtain 



R 2 (s) 



1 

2^ 



d 2 Re 



In 



(14) 



We finally carry out the differentiation and add the R\ 
contribution ([[4) to arrive at the final result 



R(s) 



1 



1 



C7T 
1 f -1 
'2^2 
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is 2 c~ 2 



(s 2 + cr 2 ) 2 (s 



-2\3 



(15) 



Eq. ( |l5| ) is applicable when the energy s » 1 and fine 
structures on scales s « 1 are inessential. Let us con- 
clude this section with a brief discussion of this result. 

The contribution R\ (cf. Eq. fll2])) has the following 
interpretation: The ^-averaged Green's function Q is 
similar to Go except for the fact that a finite width T 
has been attached to each of the uncorrelated levels £j. 
This 'smearing' implies that the corresponding correla- 
tion function Ri is Lorentzian, i.e. it is a broadened ver- 
sion of the 5-function that would be obtained for sharply 
defined autocorrelated levels. The complementary term 
i?2 describes correlations between the V degrees of free- 
dom. After combining the two contributions three qual- 
itatively different regimes can be identified: 

• For s>c _1 the dominant contribution comes from 
Ri and we obtain R(s) ~ (ttcs 2 )^ 1 . 



• c 1/2 < s < c 1 : 



Still R\ dominates but now 



Ri 



c/ir. 



3 



• s « c -1 ' 2 : The R2 contribution becomes the dom- 
inant one and we obtain the GUE-result R2 « 
— (27r 2 s 2 ) -1 corrected by a small term R\ w c/ir. 

In summary, Eq. ( |l5| ) essentially represents a superposi- 
tion of a GUE correlation function and a smeared Pois- 
sonian auto-correlation function. 



III. NON-PERTURBATIVE RESULTS 

The diagrammatic treatment is incapable of describ- 
ing structures on the energy scale of 0(A). For c> 1 an 
alternative perturbation technique, applicable over the 
whole energy axis, can be used. Within this approach 
spectral correlations are described in terms of stochas- 
tic evolution equations |p| (see also Ref. @). In this 
way one obtains a spectral correlation function R(s) that 
depends only on the combination s^/c. For large en- 
ergies, s ^> l/\/c, the result coincides with ours, i.e. 
R(s) sa (ircs 2 )^ 1 , and for small energies, s <C level 
repulsion sets in, i.e. R(s) + 1 oc cs 2 . 

The complementary regime of c < 1 can be treated 
by Efetov's non-perturbative supersymmetry technique 
p5[ (for a recent review see where averages of 

Green's functions are obtained from a generating func- 
tional. The generating functional corresponding to the 
RP model is similar to the one described in j|D • In these 
works the problem of a random banded matrix with addi- 
tional diagonal disorder was addressed. Taking the band 
width equal to the matrix size N leads to the RP model. 
From the generating functional one can obtain the cor- 
relation function of retarded and advanced Green's func- 
tions if 12 = (tr G + (E + uj)tv G~(E)) where the average 
is taken with respect to Eq. ([!]). The final integrations 
can be carried out within a saddle-point expansion the 
validity of which is controlled by A 2 //i > 1 or equiv- 
alently by c <C 1, and by w < N/fj, or equivalently by 
s <C c _1 . In the present work we skip the technical de- 
tails and concentrate on the discussion of the results. 

In the limit N — > 00 the function K 12 is given by 

K 12 (E,uj) = (tr G+(E)) (tr G~{E)) + 

+ sin (tts) + -^-j- + . (16) 

The first term is the entirely disconnected part and the 
terms of 0(c) describe deviations from a pure GUE be- 
havior. These terms represent the analogue of the contri- 
bution Ri that appeared in the diagrammatic analysis. 
They result from the correlation between the on-site en- 
ergies £j. (Note that in principle correlations of this type 
exist in the pure GUE as well. In that case, however, 
they represent negligible 0(l/A)-effect.) From ( |l6| ) we 
obtain the correlation function 




describing the spectral behavior in the regime s <C c _1 , 
c <C 1. We next turn to the discussion of this result. 
We first observe that the term c/n equals the leading 
contribution of the smeared auto-correlation i?i(s) for 
s <C c _1 . For very small level separations R(s) behaves as 
— 1 + (1 + c/n)(ns) 2 /3, i.e. apart from a slightly modified 
prefactor we obtain generic GUE-behavior. For larger 
values 1 <C s <C c~ x the leading terms are identical with 
those obtained in the diagrammatic treatment, as ex- 
pected. (By 'leading' we mean the first order terms of an 
expansion in the parameter 1/s -C 1 after the oscillatory 
structure in ( |l7j ) has been averaged out.) In particu- 
lar, the GUE behavior is only valid up to s <C l/\/c- 
Thus, the non-perturbative results underline the conclu- 
sion drawn from the diagrammatic analysis: R{s) is es- 
sentially a superposition of a GUE correlation with a 
smeared Poissonian auto-correlation. A conclusion to 
be drawn from this observation is that the analogy be- 
tween the Wigner-Dyson/Poisson transition in the RP 
model and disordered electron systems, respectively, is 
not complete. In the latter case the critical correlation 
function can hardly be interpreted as a simple super- 
position of two terms. This qualitative difference man- 
ifests e.g. in the behavior of the level compressibility, 
5 

X = lims^oo J ds R(s) (where it is essential that the 
-S 

limit N — > 00 is taken first). The two extremes GUE 
(Poisson) correspond to values x — (x = !)• In the 
case of a disordered metal at criticality the compressibil- 
ity takes an intermediate value < x < 1 01- In the 
critical RP model, however, x = 1, i.e. perturbing a 
Poisson ensemble by a GUE ensemble does not change 
the level compressibility |l3| . 

Finally, we would like to comment on the analy- 
sis [ |18|Jl9| - In these references, the RP model was solved 
exactly for arbitrary values of the parameters /1 and s. As 
a result of a sophisticated combination of supersymmetry 
and group-theoretical concepts Guhr |l8j obtained non- 
trivial double integral representations for the correlation 
functions which turned out to be difficult to evaluate. In 
order to derive closed expressions for R(s) the integral 
was analyzed in the two limiting cases c » 1 fiql and 
c <C 1 jnj by means of asymptotic expansion schemes. 
The price to pay for the mathematical rigor of Guhr's 
approach is that the physical origin of the various ingre- 
dients to R(s) is hard to identify. For this reason we 
believe, that a more conventional analysis like the one 
discussed above was calling for. 



IV. CONCLUSIONS 



We have studied the density of states correlation func- 
tion i?(s) (s measures the energy difference in units of the 
average level spacing) of the Rosenzweig-Porter model. 
This random matrix model contains a parameter /i which 



4 



allows to interpolate between GUE {n = 0) and Poisson 
statistics (jj, = oo). In the thermodynamic limit N — > oo 
the model shows three different types of universal func- 
tions R(s) depending on how fj, scales with N. From 
a diagrammatic analysis (locator expansion) assisted by 
non-perturbative methods we draw the following con- 
clusions: parameter values scaling like fi(N)/N 2 — ► 
(fi{N)/N 2 -> oo) lead to GUE (Poisson) statistics. In the 
borderline case n(N)/N 2 = c, however, a novel universal 
type of spectral behavior is observed. The corresponding 
correlation function R(s) has the following features: Like 
in the GUE-case levels repel each other, i.e. R(s) — > — 1 
for s — * 0. At some c-dependent value sq, R(s) changes 
sign, then reaches a maximum and decreases as (ttcs 2 )^ 1 
for large s. For c « 1 we find that the spectrum shows 
GUE-type statistics up to values s ~ 1/y/c. For larger 
values of s, a different type of statistics is observed, this 
being a consequence of the non GUE-correlation of the 
diagonal matrix elements in the RP model. These large 
energy correlations can be interpreted as the tails of a 
widely 'smeared' energy autocorrelation of Poissonian 
type. We thus conclude that the RP model in the critical 



case n{N)/N 2 



essentially leads to a linear superpo- 



sition of Wigner-Dyson and Poissonian behavior. Let us 
finally comment on the aspect of symmetries. In this 
paper we have considered the Rosenzweig-Porter model 
in its unitary version. It is a straightforward matter to 
extend both the diagrammatic and the non-perturbative 
analysis to the case of orthogonal respectively symplectic 
symmetry. On the other hand, none of our main conclu- 
sions on the structure of the models eigenvalue statistics 
did depend in a conceptual way on symmetry aspects. 
We thus expect the level statistics of the models of higher 
symmetry to be qualitatively similar but did not embark 
on any kind of quantitative analysis. 
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Figure 1: (a) Graphical representation of the Green's 
function Go, the perturbation matrix elements Vij and its 
correlations. (b),(c) The Dyson equation for the full and 
the averaged Green's function, respectively, (d) Graphi- 
cal representation of the self-energy, (e) A diagram with 
crossed GUE lines, (f) Explanation see text. 

Figure 2: Graphical representation of S n (E,E') appear- 
inginEq. @. 
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